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A method is proposed for manipulating with diagrammatic expansion of 
Green's function of Frenkel's exciton random walks on the perfect lattice. The 
method allows one to select diagrams, to supply diagrams with factors con- 
taining information about the number of sites the diagram has passed through, 
etc. Simple problems related to the defect lattices are considered using the 
proposed method. The new criterion of localization of Frenkel exciton - the 
number of sites covered by the wave function - is established. The number of 
sites covered by the zero state of ID non-diagonally disordered chain is stud- 
ied. It is shown that this problem can be solved by calculating the random 
walks Green's function with modified diagrammatic expansion. By means of 
the developed method an exact analitical expression for the number of sites 
covered by the zero-state is obtained and zero-state is shown to be localized. 



Let us consider some lattice with sites ri...r]\f and introduce matrix W with elements 



The function w{r) is supposed known. The problem of Frenkel exciton random walks 
along the lattice ri...rN of twolevel atoms, can be reduced to the matrix of this kind [0 and 
this problem is not the only one. In this case w{r — r') represents an interaction which is 
able to transport an excitation from atom r to atom r' . The study of Green's function (GF) 



I. INTRODUCTION 
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which corresponds to the above matrix (1) is of particular interest We write down this 
function in convenient dimensionless form: 

G = (l-W)-i (2) 

If suppose the summation over the repeating indexes, the following series for GF may be 
written: 

'^l"^ = (t—^) = ^Irn + Wlm + WlrWrm + WwWrpWpm + ••■ (3) 
V i — \\ J Im 

If we denote matrix element W^p by an arrow directed from site r to site p , then each 
product of matrix elements in (3) is being represented by diagram (trajectory) , which 
connects sites / and m . Then (3) can be rewritten as: 

Gira = ^ Sum of all diagrams, which connect sites / and m ^ (4) 

The similar method has been used in 0. The above representation allows one to write 
down the formal expressions for Green's functions of a certain relative problems. For ex- 
ample, if certain atoms ri...rt have been removed, then the expression for Green's function 
must be modified as follows: 

(Sum of all diagrams, which connect sites / and m, \ 
(5) 
not passing through the sites ri...rt J 

If we now deal with lattice whose sites can be occupied by atoms with probability C < 1 
and we are interested in the averaged Green's function, then each diagram D of this function 
should be multiplied by probability that all n{D) sites D has passed through, are occupied 
by atoms. This probability is equal to C"^^^ and, consequently, we can write for the averaged 
Green's function : 

/ Sum of all diagrams D, which connect sites / m, \ 



{Glm)c 



(6) 



each diagram is being multiplied by C"'-^-', 
\ where n{D) - is the number of sites in the diagram J 
In the end if we are interested in {Gim)c, C ~ 1, we have to calculate d{Gim)c/dG , 
C = 1. Making use of (6), it is easy to obtain: 



/ Sum of all diagrams D, which connect sites I and m, \ 



d{Gim)c 

dC 



Each diagram is being multiplied by n{D), 
where n{D) - is the number of sites in D 



(7) 



V 



These examples shows that Green's functions, related to various problems, can be ob- 
tained by manipulating with one and the same series of diagrams. So, it is desirable to 
construct some kind of observer who, while mov ing along the diagram, would register the 
sites which diagram has visited, would count the number of sites which diagram has passed 
through and so on. For example, in order to obtain Green's function (5) the above observer's 
duty should be as follows: While moving along the diagram and while visiting one after an- 
other the sites of this diagram, observer should check whether site he just has visited is one 
of ri...rt which have been removed and, if it is so, observer should multiply diagram by zero. 

To obtain Green's functions (7) and (6), observer should count sites which he has passed 
through and should multiply the diagram by total number of sites n{D) in the diagram or 
by C"^^^ respectively. In section 2 we show that in certain cases it is possible to build up the 
mathematical realisation of described observer by replacing matrix elements Wrp in Green's 
function (4) by specially constructed watching operators. 

However, the random walks problems described in section 2 which can be solved up to 
the very end (at least at present) by above method is rather simple and are of some interest 
only as a demonstration of absence of errors in the mathematical scheme of the method. 

In section 3 we consider much more nontrivial problem related to the character (localized- 
delocalized) of Frenkel's exciton zero-state in ID disordered chain. We establish a novel 
criterion of localization - the inverse number of sites covered by the wave function - and 
show that the problem of calculation of this value may be formulated in terms of diagrams 
and exactly solved by means of specially built watching operators. 
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A. 



Let us consider the Green's function Gor which is defined on certain periodical simple 
lattice. Let us now perform the Fourier transformation with respect to its second index 

9k = ^ r , where fk = V w(r) exp(zfcr) (8) 

^- Jk r 

in terms of diagrams: 

/ Sum of all diagrams which are starting from the site 0. \ 
gk — Each diagram is being multiplied by exp{ikrf), where r/ - (9) 
\ is a coordinate of the last site in the diagram / 

Let us consider formula (5) and denote by L the set of sites occupied by atoms and by 
T the set of sites ri...rt. 

In order to convert function (9) to the Green's function of this unperfect (defect) sys- 
tem we replace matrix elements Wrp in (9) (which are denoted by diagram's arrows ) by 
the operators Wrp which act on upon arbitrary function Jq {q is an auxiliary variable) in 
accordance with the rule: 

WrpJq = E - qy'^'~''^J,': where K{q) ^ ^ ^M^QR) (10) 

Summation over q is carried out over the first Brillouin zone. Now Green's function (9) 
becomes the operator gi- — > g;.. Let's apply this operator to unit. 

Let's consider a certain diagram and apply the operator which corresponds to the first 
arrow to unit- Worl- It is easy to see that if r e L, then the result will be H^or exp(zgr). 

Otherwise we obtain zero . Consequently, if the first arrow of the diagram came to the 
empty site, the diagram is being multiplied by zero. By induction we come to the conclusion 
that it is a general result - if certain arrow of the diagram reaches the empty site p (i.e. p 
is one of the ri...rt) the operator corresponding to this arrow gives zero. 

We see that operators (10) give the mathematical realisation of the observer described 
in the end of the previous section. That's why we will call them the watching operators. 



dk+q (11) 



Consequently, if we denote the Green's function corresponding to the lattice with defects by 
dg , the following formula may be written: 

/ Sum of all diagrams which are started from site 0, \ 
which are not passed through ri-.-r^. 
Each diagram is multiplied by exp{i{k + q)rf), where r/ 
\ - the coordinate of the last site of the diagram / 
By replacing Wor — w{r) in the explicit expression (8) for gk (9) by the watching operators 
(10) and by setting A; = 0, we obtain the equation for dg-. 

1 = (1 - fo)dg, where /o = ^ Wor (12) 

r 

The simple calculations shows that dg may be introduced in the form: 

J2 exp{iqR) 
dq = ""^^^ _ , (13) 

-•- Jq 

where parameters Xr may be obtained from the equations: 

E Gr-r^Xr, = Gr, Re T, (14) 
R'eT 

where 

GR^^j:e-"'V{i-f,) (15) 

q 

- is the Green's function of the periodical lattice. If the amount of removed atoms is not 
too large, the solution of (14) may be presented in a compact form . For example, if the only 
one atom have been removed from the site R, then the solution of (14) is: Xji = Gr/Gq. 

The further analysis of (14) is possible either if the empty sites forms the periodical 
structure. It should be pointed out that while obtaining (14) we have supposed that the 
site G L. If it is not so, then dg = 0. Formulas (13) and (14) gives it (in this case 
Xq — l,Xji^Q — 0), but an exact solution of (12) is non-zero and represents the sum of 
diagrams which start from the zero site and never return back to this site more. 

At last let us note that the initial site is not especial one and the Green's functions 
having another first index may be calculated by the same way. 



B. 



In this subsection we will introduce another type of watching operators which give possi- 
bility to multiply the diagrams by factors, containing the information about the sites which 
diagram has visited. 

Let's introduce the watching operators dependent on a parameter a which operate as 
follows: 

WrpJ, ^ Wr,(l + ^E)^''^"'"^-^^ (16) 

Let's apply the operator corresponding to certain diagram D 

containing in to some function $g = X]ii0i?exp(— zgi?) (summation over all lattice). 
Operator corresponding to the first arrow of the diagram operate as follows: Wor-^g — 

Wor exp(2gr)$^, where = J2R<^'RCxp{—iqR), and 0'^ = (pR\i R ^ r and 0^ = 0:0^. 

Moving further along the diagram and applying step by step the watching opera- 
tors corresponding to the arrows of the diagram , one can see that the diagram D will 
be multiplied by exp{iqrf{D))^g{D) = exp{irf{D))J2R'pRiD)cxp{—iqR), and if i? 7^ 0, 
MD) = (pRa^^^"" of of t^^^ ^if^ ^ by diagram D)_ ^^^^ _ ^^^^ ^-^^ 

the diagram D. 

The site ( which the diagram started from) is being taken into account in a special 
way: 0o(^) = ^oa^™^''^ of diagram's D visits of site -1) Consequently, %{D) carry the 
information about sites, which the diagram D has visited and about corresponding number 
of visits. For this reason we call this function carrier. The arbitrary function which we 
have acted on upon by operator, may be called the carrier in the initial state. So, we see, 
that after applying gk (containing the watching operators (16)) to the carrier in initial state 

each diagram is being multiplied by her own carrier. Let's now obtain the function (7) 
by means of the watching operators (16). It is clear now, that if we take the carrier in initial 
state in the form: $q = J2rj^o exp(— zgi?) + a — NSgo + a — 1 and apply gk (containing the 
watching operators (16)), the following formula is justified: 
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9k^q 



7. (17) 



/ Sum of all diagrams which are started from site 0. \ 
Each diagram is being multiplied by carrier: X)/? 6~*''^q;^*^'^'^\ 
Where V{D, R) - the number of diagram's D visits of site R. 
Each diagram is being multiplied by e^^'^'^'^^'^f , 
\ where Tf - is the coordinate of the last site of the diagram. J 

The equation for 7^: $q = (1 — fk)'yq ( being built by using the watching operators (16)) 
may be easily solved in terms of Green's functions (8) and (15): 

(« - '^)9q+k9k 



Iq = NgkSgo + 



(18) 



l-(«-l)(Go-l) 

The comparison of the diagrammatic formulas (7) and (17) shows, that the Fourier 
transformation of (7) with respect to the second index (similar to that in (8)) may be 
represented in the form: 

d{gk)c 



(19) 



An other types of the watching operators are also possible @]. In the next section we 
apply the developed method for solving the problem which considerably differs from those, 
which we have considered in this section. 



III. THE EXPLORATION OF FRENKEL EXCITON ZERO-STATE BY MEANS 
OF THE WATCHING OPERATORS METHOD 

A. The criterion of localisation 

Let us consider the random process on the ID chain of sites i = 1,2, 3... n, which is 
defined by the following recurrent relation: 

^Pi = r]iLPi_i, ipo = I (20) 

Where rji = Ui/U-, and Ui , Ul - are independent, similarly distributed, positive random 
variables, which can have the values a > and 6 > with probabilities C and 1 — C 
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respectively {a > b). (pi represents the square module of the Frenkel's exciton zero-state in 
the ID non-diagonal disordered chain 0, when the Hamiltonian matrix elements can have 
values y/a and y/b with probabilities C and 1 — C respectively. Let us characterize the zero- 
state by the number of sites A^*, covered by the corresponding wave function and calculate 
this number as follows. If = on certain site i , then this site is not covered and is 
being not taken into account while calculating N*. Vice versa, if (pi reaches its maximum 
ipmax on the site i, then this site is covered completely and its contribution while calculating 
A^* is equal to unit. In the general case, the contribution of the arbitrary site i is equal to 
^i/^max and we obtain the following formula: 

N* = (21) 

may be treated as a criteria of localisation - it is finite when zero-state is localized 
and is going down to zero as the inverse chain length n~^, if zero-state is delocalized. 

If lim A^* < oo, then for the averaged quantities the following formula may be written: 

n— >oo 

lim (A^*) = ((A^*)) (22) 

Where () denotes the averaging over the realisations of the random process and (()) 
denotes the averaging over all chains lengths n. In the next subsection we will describe the 
method for exact calculation of {{N*)) by means of methods developed in section 2. 

B. The diagrammatic representation of the random process ipi. The method for 

calculating {{N*)) 

Let us introduce the diagrammatic representation of random process (pi as follows. The 
four possible values are acceptable for ipi. Each value may be represented by horisontal and 
vertical arrows as follows: 



^In the considered random process always (pi ^ 0, but it is not important for this reasoning 
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' a/a t 
h/h i 
a/h - 

Ma - 



CH 



{I -eft 

-- (1 - C)Ct 
-- (1 - C)Ct 



(23) 



The length of each arrow is supposed to be equal to unit. Then the arbitrary reahsation 
of random process <^j, related to the chain consisting of n sites, is being represented by a 
trajectory - diagram - on the two dimensional square lattice. These diagrams has n arrows. 
Without loss of generahty we can suppose that each diagram is being started from the zero 
site of this lattice. 

Now let us ascribe to the i-th arrow of the diagram the probability of the corresponding 
value of ?7j, in the same manner it have been done in (23). The parameter i , < i < 1 -is 
an auxiliary one and in the end of the calculation t — > 1. 

Consequently, the diagram completely describes the distribution of r]i in corresponding 
realisation and the diagram's value is equal to the probability of this realisation. Let's now 
explain how the (pmax and can be extracted from the diagram. 

Introduce the coordinate system XY with zero point coinciding with the zero of the 
above mentioned square lattice. 

If r is the x- projection of the right extremal site (sites) of the diagram, then ipmax = C^, 
where ^ = a/b > 1. 

In order to calculate J2i for certain realisation, let us suppose that we have multiplied 
the corresponding diagram D by the following carrier ^q{D) = J2x ^(-D, x) exp{—tqx) , where 
V{D,x) - the number of diagram's visits of the sites whose x- projection (contrary to (17)) 
is equal to x. 

After the Fourier transformation the carrier becomes: (j)x{D) = V{D,x). Then J2i^i = 



Now let's introduce the matrix W (1) defined on the above mentioned square lattice by 
means of the following function w{r): 
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w{r) — < 



(24) 



' CH, if r = (0, 1) t 
(l-C)H, if r= (0,-1) i 
(l-C)Ct, if r= (1,0) ^ 
[(l-C)Ct, ifr=(-l,0) ^ 
and build up the following function in accordance with the first section: 

f Sum of all diagrams which are starting from the site 0. \ 
Each diagram is being multiplied by the following carrier: 

where V{D,x) - is a number of diagram's D visits 
of the sites whose x- coordinate is x. 
The diagram is also multiphed by 
^i{k+q+q')xf ^ where Xf - x- coordinate 
\ of the last site of the diagram / 

Defining the operation D/Dx as Dy/Dx = y{x + 1) — y{x), we introduce the function: 



F(q,q',k,a,(3) 



(25) 



Where iV - is the number of sites along the x- direction. — cxd in the end of the 
calculations. In terms of diagrams the expression for U {x, x') is as follows: 

/ Sum of all diagrams which are being started from \ 
the site 0. All diagrams are being multiplied by the carrier 
{5x,r-5x,i-^V{D,x'), 
U{x, x') — where V{D, x') - is the number of diagram's D visits (27) 

of the sites whose x- coordinate is x. 
Where r{l) - x- projection of the right (left) extremal site 
\ (or sites) of the diagram D. J 

Let us explain the sense of operation D/Dx in (26). After the Fourier transformation 
with respect to q and integration over a, the dependence on x for carrier of any diagram is 
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being represented by a set of ordered S- symbols: 5x,r + Sx,r-i + Sx,r-2 + ■■■5x,h where r{l) - 
is the X- projection of the right (left) extremal site (or sites) of the diagram. 

The absence of gaps in the series of the S- symbols is take place because an arrows, 
which correspond to the function ^(r) (24), can connect only neighbour sites (along the x- 
direction) 

Consequently, after operation D/Dx only right 5x,r and left 5x,i-i 5- symbols (shifted 
per unit) are survive. Then the following formula for the average number of sites (22) is 
justified: 



oo oo 



((7V*))=lim(l-t);^ ^ U{x,x'r'-^ (28) 

X = x' = — OD 

The 1 — t - is the normalisation factor related to the averaging over chain lengths and 
is equal to inverse sum of all diagrams, i.e. g^^ (formula (8), constructed by means of 
(24)). The exphcit expression for the function (25) can be obtained by manipulating with 
diagrammatic expansion of the Green's functions (8), (9), by means of the watching operators 
of the following form: 



rr=(0,l) C^t[l + s^V ][i + ^E,/ 



^(r) 



r=(0,-l) (1-C)2t[l + ^EJ[1 + ^E,'] 

r = (1, 0) (1 - C)Ct[l + ^ ] [1 + y E,' ]e<''+^') 



(29) 



'=(-1,0) (1 - C)Ct[l + ] [1 + y E,' ]e-^(^+«') 

These operators differ from (16) by that depend on two auxiliary wave vectors g, q'{—Tr < 

q,q' < tt) , and that they watch only x projections of the diagrams arrows on two dimensional 

lattice. By applying to {NSgo + a — l){N5q'Q + /3 — 1), we obtain the standart equation 

for F{q, q', k, a, (3): 

(1 - fk)F = {N5,o + a- l){N5,,o + /3 - 1) (30) 

The solution of this equation is rather clumsy. It may be reduced to the integral equation 
of convolution type, which can be solved by Fourier transformation. The solution can be 
expressed in terms of known Green's function of the Frenkel's exciton on one dimensional 
lattice with only nearest neighbour interaction. 
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Using (25), (26), (28) we obtain the exact expression for A^*: 

Let us to note in the conclusion that the above problem is not the most simple problem, 
related to the random process (20), which can be solved by means of the watching operators 
method . By applying the diagrammatic representation of the random process (20) and 
the watching operators method one would obtain such nontrivial characteristics of (20) as 
distribution function for if) max and the dependence of {ipmax) as a function of chain's length 
n. 

We present the result (31) as an argument in favour of localized type of the zero state, 
contrary to ||^, where zero state was declared to be delocalized because of its zero inverse 
localization length (ILL). This criterion of localization is not correct because the above 
model also has zero ILL, but (31) shows that zero state is always localized. 



[1] LM.Lifshits, S.A.Gredeskul and L.A.Pastur, Introduction in Theory of Disordered Systems, 
Nauka, Moskow (1982). 

[2] C.R.Gochanour, H.C.Anderson, M.D.Fayer, J. Chem. Phys., 70 (9), (1979) 4254. 

[3] G.Theodorou, M.H.Cohen, Phys. Rev. B 13 (1976) 4597. 

[4] Kozlov G.G. " On the Green's function of disordered sysytems" , presented to VINITI, 1993, N 
DO-85-79. 



12 



